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Abstract 

A direct approach to exact solutions of nonlinear partial differential equations is 
proposed, by using rational function transformations. The new method provides a more 
systematical and convenient handling of the solution process of nonlinear equations, 
unifying the tanh-function type methods, the homogeneous balance method, the exp- 
function method, the mapping method, and the F-expansion type methods. Its key point 
is to search for rational solutions to variable-coefficient ordinary differential equations 
transformed from given partial differential equations. As an application, the construction 
problem of exact solutions to the 3 + 1 dimensional Jimbo-Miwa equation is treated, 
together with a Backlund transformation. 
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1 Introduction 

Partial differential equations (PDEs) describe various nonlinear phenomena in natural and 
applied sciences such as fluid dynamics, plasma physics, solid state physics, optical fibers, 
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acoustics, mechanics, biology and mathematical finance. Their solution spaces are infinite- 
dimensional and contain diverse solution structures. It is of significant importance to solve 
nonlinear PDEs from both theoretical and practical points of view. Due to the nonlinearity 
of differential equations and the high dimension of space variables, it is a difficult job for us 
to determine whatever exact solutions to nonlinear PDEs. 

It has been a successful idea to generate exact solutions of nonlinear wave equations by 
reducing PDEs into ordinary differential equations (ODEs). Many approaches to exact solu- 
tions in the literature follow such an idea, which contain the tanh-function method the 
sech-function method [6j-[9], the homogeneous balance method fiUl [TT] the extended tanh- 
function method y^-^Mj, the sine-cosine method [T3], the tanh-coth method [TU], the Jacobi 
elliptic function method [TTj, the exp-function method [H], the F-expansion method [19j, 
the mapping method [20], and the extended F-expansion method [211 [22] • Given an ODE of 
differential polynomial type, either constant-coefficient or variable-coefficient, one can always 
adopt computer algebra systems to search for rational solutions pretty systematically. This is 
one of the main reasons why those reduction methods work well. 

Based on this observation, we would, in this paper, hke to propose a direct and systematical 
approach to exact solutions of nonlinear equations by using rational function transformations. 
The method is very suitable for an easier and more effective handling of the solution process 
of nonlinear equations, unifying the existing solution methods mentioned above. Its key point 
is to find rational solutions to variable-coefficient ODEs transformed from given nonlinear 
PDEs. Together with an auto-Backlund transformation, an application of our approach will 
generate new exact solutions to the 3 + 1 dimensional Jimbo-Miwa equation [23] (called the 
Jimbo-Miwa equation in [21]): 

-P,/m(^) •= Uxxxy + 3liyZ/a;a; + SUxU^y + 2Uyt — Su^z = (1-1) 

where u = u{x,y, z,t), Ux = etc. 

The 3 + 1 dimensional Jimbo-Miwa equation fll.ip is the second member in the entire 
Kadomtsev-Petviashvili hierarchy [23], originally defined by a Hirota bilinear equation 

[{Dl + 2Dt)Dy - 3DxDz)]r.r = 0, (1.2) 

with the link being u = 2(lnr)2;. The Jimbo-Miwa equation f 1 1.1 1) passes the Painleve test 
only for a subclass of solutions [21] and its symmetry algebra does not have a Kac-Moody- 
Virasoro structure [2E|- Nevertheless, different type solutions to the Jimbo-Miwa equation 
(11. ip are found (see, say, [26]- [30]). The Hirota perturbation technique yields one- and two- 
soliton solutions [21] and dromion type solutions [2H] , and the truncated Painleve series leads 
to a special Backlund transformation [31]. Obviously, the Jimbo-Miwa equation (II. ip has the 
following X- or y-independent solutions: 

u = f{x, t) + h{z, t), u = g{y, z) + h{z, t), (1.3) 

where /, g and h can be arbitrary functions in the indicated variables. These solutions contain 
more special solutions: u = f{x,t), u = g{y,z) and u = h{z,t), which are independent of 



2 



two variables of x,y,z,t. Starting with such solutions, various variable separated solutions 
have been presented by using the truncated Painleve series [HT] [32] and abundant nonlinear 
coherent structures have been exhibited [32]. We will concentrate primarily on constructing 
travelling wave solutions to the Jimbo-Miwa equation f 1 1.11) by formulating a problem of finding 
rational solutions to a transformed Jimbo-Miwa equation. 

The paper is organized as follows. In Section [2], a unified formulation for getting exact 
solutions to nonlinear equations is proposed, by using rational function transformations. In 
Section [3l an application and an auto-Backlund transformation are made to solve the 3 + 1 
dimensional Jimbo-Miwa equation (II. ip . A few of concluding remarks are given in the final 
section, along with some polynomial solutions. 

2 A transformed rational function method 

To describe our solution process, let us focus on a scalar 3 + 1 dimensional partial differential 
equation 

P{x, y, z, t, u^, Uy, Uz, ut,-- ■) =0, (2.1) 

though the solution process also works for systems of nonlinear equations. We assume that 
there are exact solutions to the differential equation fl2.ip : 

u{x, y, z, t) = ^ = ^{x, y, z, t). (2.2) 

Usually, we can have 

^(x, y, z, t) = ax + by + cz — ut, (2.3) 
where a, b, c and u are constants, in the constant-coefficient case, and 

C{x, y, z, t) = a{t)x + b{t)y + c{t)z - uj{t), (2.4) 

where a(t), b{t), c{t) and uj{t) are functions of t, in the t-dependent-coefficient case. Under the 
transformation (12.21) . the partial differential equation (12.11) is put into an ordinary differential 
equation: 

Q{x,y,z,t,u^^\u^^+'\---)=0, (2.5) 

where = ^ > 1, and r is the least order of derivatives in the equation. To keep 
the solution process as simple as possible, the function Q should not be a total ^-derivative 
of another function. Otherwise, taking integration with respect to ^ further reduces the 
transformed equation. 

An important step in the solution process is to introduce a new variable r] = rjl^) by a 
solvable ordinary differential equation, for example, a first-order differential equation: 

rj' = T = Ti^,r)), (2.6) 

where T is a function of ^ and rj, and the prime denotes the derivative with respect to ^. In 
case that we have a general second-order differential equation to begin with, we should first 
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obtain its first integrals [33] and then use the method of planar dynamical systems to solve 
[51] . Two simple solvable cases of the above function T are as follows: 



T = Tij]) = rj, T = T{r]) = a + •rf , a = const. (2.7) 
The corresponding first-order equations have the particular solutions r] = and 



V 



— |, when a = 0, 

— \/—a tanh-\/— or — \/ —a coth \/ —a ^, when a < 0, (2.^ 
^/a tan i/a or — ^/a cot -^/a when a > 0, 



respectively Those two cases correspond to the exp-function method and the extended 
tanh-function method, respectively. 

More general assumptions than fl2.7p can engender special function solutions to nonlinear 
wave equations. For instance, taking (77')^ = Sir]) with some fourth-order polynomials Sir]) in 
?7 (or equivalently, 77" = Rij]) with some third-order polynomials Rij]) in 77) can yield Jacobi 
elliptic function solutions; and such assumptions are the bases for the extended tanh-function 
method, the F-expansion method and the extended F-expansion method, and work for many 
particular nonlinear wave equations. 

The basic idea of using solvable ordinary differential equations was successfully used to solve 
the 2+1 dimensional Korteweg-de Vries-Burgers equation in [33], based on ar]"+br]'+c'r]'^+dr] = 
{a,b,c,d = const.), and the Kolmogorov-Petrovskii-Piskunov equation in [T2], based on 
f]' = 1± rf. Later it was broadly adopted in the extended tanh-function method [131 115 5 the 
tanh-coth method [16], the F-expansion method [I9], the mapping method [201 ES], and the 
extended F-expansion method [211 [22] • 

Let us proceed to consider rational functions 

where m and n are two natural numbers, and pi, < i < m and qi, < i < n are nor- 
mally constants but could be functions of the independent variables as in the F-expansion 
type method. All Laurent polynomial and polynomial functions are only special examples of 
rational functions. We search for travelling wave solutions determined by 

n«(0=^(r/) = ^, (2.10) 

where pi^r]) and qi^r]) are polynomials as indicated above. It is direct to compute that 



u 



(''+^) - Tv', = Tdr,u^''^ = Tv" + T'v', ■■■ , (2.11) 



which is based on = Tdrj. Note that by the prime, we denote the derivatives with respect 
to the involved variable, for instance, u' = v' = and v" = 
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Now we assume that the transformed equation (12.51) is a rational function equation of r] 
with a given pair of m and n. This can be achieved for all nonlinear equations of differential 
polynomial type, when T is, for example, a rational function in rj. Then we just need to force 
the numerator of the resulting rational function in the transformed equation to be zero. This 
yields a system of algebraic equations on all variables a, b, c,u, pi, < i < m and Qi, < i < n, 
and solve this system (which could be a differential system as in the F-expansion type method) 
to determine pi^r]), qiji) and ^. Finally, integrating virf) with respect to ^, r times, we obtain 
a class of travelling wave solutions 

u{x,y,z,t) = u{0= [ ■■■ [ ^W^dC---dC 



where di, 1 < i < r, are arbitrary constants. If r = 1, there is only the last definite integral 
over [0,^] in (I2.12p . The resulting solutions will definitely contain a polynomial part in ^, 
when r > 1. 

The case of the exp-function method: 

If we take r] = e^, then the solution function is determined by 

g^e"? + qn-ie^"" ^'^ ^ h go 



This can yield solutions generated by the exp-function method 
The case of the extended tanh-function method: 

If we take i] = — tanh((^), then the solution function is determined by 



I I (-l)>^tanh-(0 + (-1)— Wi tanh— ^(0 + ■ ■ ■ + po 
"^^^ ' -l)"g„tanh"(0 + (-l)"-ig„-itanh"-^(0 + --- + go 



This can yield solitary wave solutions. If we take 77 = tan(^), then the solution function is 
determined by 

u(n = I ! P"^ tan"-(0 + Pm-i tan"--i(0 + ■ ■ ■ + Po 
' ' gntan"(0 + gn-itan-i(^) + --- + (/o 



This can yield periodic wave solutions. The other choices of 77 = — coth(^) and r] = —cot{^) 
generate similar type exact solutions. Note that we only use the selection of a = 1 and 
a general nonzero value of a may lead to more general solutions. The selection of a = 
presents rational solutions. If v is a polynomial and r = 0, then u gives an exact solution in 
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the form of a finite series in tanh,^ or tan^, which is obtainable by the tanh-function type 
method PIHIIIS]. 

The solution process described above unifies the existing methods using tanh-function 
type functions, tan-function type functions, the exponential functions and the Jacobi elliptic 
functions, and allows us to carry out the involved computation more systematically and con- 
veniently by powerful computer algebra systems such as Maple, Mathematica, MuPAD and 
Matlab. While applying to construction of special function solutions to nonlinear equations, 
we have to pay special attention to the particular forms of given nonlinear equations to get 
workable transformed equations. 

3 Solving the 3+1 dimensional Jimbo-Miwa equation 

To generate travelling wave solutions 

u{x, y, z, t) = u{^), ^ = ax + by + cz — ut, (3.1) 

where a, b, c, are the angular wave numbers and u are the wave frequency, we only need to 
solve the reduced 3 + 1 dimensional Jimbo-Miwa equation f 1 1.1 1) : 

a^bu^^^ + Ga^bu'u" - {2bu + ?,ac)u" = 0, (3.2) 

where the prime denotes the derivatives with respect to ^. Integrate it once with respect to ^ 
to obtain 

aV + 3a^b{u'f - {2bu + 3ac)u' = 0. 
We set r = 1 and u' = v, and then, we have the transformed Jimbo-Miwa equation 

a^T^v" + a^bTT'v' + 3aV - {2bu + 3ac)v = 0, (3.3) 

where the prime denotes the derivatives with respect to rj. 

3.1 The case of r]' = r] 

In this case, the transformed Jimbo-Miwa equation (13. 3p becomes 

a^br]'^v" + a^br]v' + ?,aV - {2bu + 3ac)v = 0. (3.4) 

We try to search for a rational solution v with m = n = 3. A direct computation with 
Maple tells that there are only two choices with a non-constant v and abc ^ 0: 



4:aqiq2'i] a{a^b — 3c) 

Aq^rf + Aqiq2ri + ql' 2b 
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and 

, . _ PojpW + l6poPiV + l6pl) _ 3po _ 9po(3bpo + egg) 
qoiPlV - 8poPiV + ^QPo) Qo 26gg 

Accordingly, we have the following travelling wave solutions to the 3 + 1 dimensional Jimbo- 
Miwa equation (13.31) : 

u{x, y, z, t) = — - — - + d, ^ = ax + by + cz — 1, (3.7) 

qi + 2^26? 2b 



and 



?/, ^, ^) = ^1 TT + — C + (i, ^ = X + by + cz — ^ 1, {3.i 



where the involved constants are all arbitrary. Note that the second class of travelling wave 
solutions, defined by (13.81) . contain a linear function in ^. 

We point out that all exact solutions generated in [29] belong to the above first class of 
travelling wave solutions, defined by (13.71) . which can be observed by multiplying common 
factors, re-scaling the involved constants and checking the rational form of the ti-function. 
Only for the solution in the formula (30) of [22!, ^^^^ needs to cancel one common factor 
e« + bi. 



3.2 The case rj' = a + rf 

In this case, the transformed Jimbo-Miwa equation (13. 3p becomes 

a^hq'^v" + 2a%rfv' + 2a^barfv" + a^ba^v" + 2a^bar]v' + 3a^bv^ - {2bu + 3ac)v = 0. (3.9) 

We try to search for a rational solution v with m = 3 and n = 1. Similarly, a direct 

computation with Maple tells that there are only two choices with a non-constant v and 
abc 7^ 0: 

, , 9 aiAba^a + 3c) , , 

vM = -2aa - 2arf , u = ^ ; 3.10 

2b 

and 

, , 2 9 a{4:0?ba — 3c) , , 

v{n) = --aa - 2arf, u = ^ '-. (3.11) 

Through re-scaling the involved constants, we can know that the general nonzero value of 
a does not engender more general solutions in the above two cases. Therefore, we select a = 1 
and take r] = tan^ and 1] = — cot^, we obtain the corresponding travelling wave solutions to 
the 3 + 1 dimensional Jimbo-Miwa equation (13.31) : 

/ N ^ r ^ , aUba^ + 3c) 
u{x,y, z, t) = —2a tan 4 + a, t, = ax + by + cz -\ — 1, 

. (3.12) 

4a^ a{4a%-3c) 
u{x, y, z, t) = —2a tan t, H — - — ha, ^ = ax + by + cz — 1, 
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and 

/ \ ^ , ^ , a(4ba'^ + 3c) 
u{x, y, z, t) = 2a cot c, + a, t, = ax + by + cz -\ 1, 

, (3.13) 

X . . 4a^ , . , a(4a%-3c) 

uix, y, z, tj = la cot t, H h ct, t, = ax + by + cz ; 1, 

3 2b 

where the involved constants are all arbitrary. Note that each second class of travelling wave 
solutions contain a linear function in ^. The selection of a = leads to a class of rational 
solutions 

2a Aab 
y^^^') = J + d= 2abx + 2b-y + 2bcz + 3act + ^' ^'-^'^ 
where the involved constants are all arbitrary. 

It is also direct to see that if we take rj = — tanh ^ and rj = — coth ^, we will obtain the same 
travelling wave solutions as those presented in the sub-section 13. 1[ Our solutions generated 
from (13. 5p and f l3.10p contain the solutions given in [27j, and reduce to the solutions in the 
2 + 1 dimensional case oi u = u{x, y, t) in 



3.3 Backlund transformation 

Let u = u{x,y,z,t) be a solution to the 3 + 1 dimensional Jimbo-Miwa equation (11.11) . Evi- 
dently, if a function v = v{x, y, z, t) satisfies 

3uyVxx + ^UxxVy + Su^yVx + 3uj:Vxy + Pjm{v) = 0, (3.15) 

where Pjm is the Jimbo-Miwa operator defined in (11. ip . then the sum of the two functions, 
w = u + V, gives another solution to the Jimbo-Miwa equation (II. ip . Therefore, once we find 
a function v satisfying (I3.15p . we get a new solution w = u + v from a known one u. This 
forms a general auto-Backlund transformation for the Jimbo-Miwa equation (11.11) . 

Sub-BTl: First, it follows directly from the above Backlund transformation that if two 
solutions u and v of the Jimbo-Miwa equation (ll.ip satisfy 

'^yVxx ~\~ IJ'XxVy + UxyVx + UxVxy = 0, (3.16) 

then w = u + V is a, third solution to the Jimbo-Miwa equation (II. ip . Further, we easily see 
that any solution u = u{x,y,z,t) plus an arbitrary function h{z,t) gives a new solution for 
the Jimbo-Miwa equation (II. ip . 

Sub-BT2: Second, if we take 

^; = 2(ln0),. = 2^, (3.17) 
then w = u + 2(ln (j))x solves the Jimbo-Miwa equation (II. ip . when (p = (p{x, y, z, t) satisfies 

'Pxxxy + "^Ux^txy + 'iUy't'xx + 20„f — "icfixz = 0, 

(3.18) 

(Pxxx(py + ^(Px(pxxy " ^(pxx(pxy + 'iUx(j)x4>y + 3ttj^(;/)^ + 20^0^ - ?>(t)x4>z = 0, 
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where w is a solution to the Jimbo-Miwa equation (11. ip . This special Backlund transformation 
was introduced through the truncated Painleve series for the Jimbo-Miwa equation (11.11) in [31] 
and for the derivative Jimbo-Miwa equation {Pjm{u))x = in [38]. It was successfully used 
to present diverse variable separated solutions involving arbitrary functions of two variables 
for the Jimbo-Miwa equation in [31] |32j . 

One one hand, if we take a solution u = f{x, t) + h{z, t) and 

3 3 
= ao + aip{—y + z) + a2q{x, t) + asp{—y + z)q{x, t), (3.19) 

Co Co 

where a,, < i < 3, and co are constants, then the two conditions in (I3.18P will be satisfied 
when 

2gt + qxxx - qx{co - 3/^;) = 0, (3.20) 

which determines the function /. All the involved constants and the functions h,p,q are 
arbitrary. This generates the class of exact solutions to the Jimbo-Miwa equation (II. ip . 
presented in [32]. More specially, when ai = 02 = 0, we obtain a sub-class of solutions given 
in 



On the other hand, if we take a solution u = g{y, z) + h{z, t) and 

(p = ai + a2p{y,z)q{x,t), (3.21) 
where ai and 02 are constants, then the two conditions in (13.180 will be satisfied when [31] 

{qxf - qqxX = 0, Pyiqxxx + '^qt) + Sgypqxx - ^Pzqx = 0. (3.22) 
A particular case is given by 

q{x, t) = e-+-*, ,(,, z) = -^ [ + ^'^^Py - ^'^P^ dy, (3.23) 

J p 

where Ci and C2 are constants, and thus the resulting new solution reads 
2a2Cip(i/,2)e^i^+^^* 1 f (4 + 2c2)py - 3cip^ 
" = a, + a2piy,z)e^^^^^^^ ' M J ~P + ^' ^ 



where all the involved constants and functions are all arbitrary. 

Sub-BT3: Third, if we take a travelling wave solution u = u{x, y, z, t) = u{ax+by+cz—u}t) 
to the Jimbo-Miwa equation (11.10 . then the function 

w{x, y, z, t) = u{ax + by + cz — ut) + a'x + h'y + c z — to't + d', (3.25) 

where a', b', c', d', u' are all constants satisfying 

ab' + a'b = 0, (3.26) 

presents a new solution to the Jimbo-Miwa equation (11.10 . 

Now taking advantage of this feature in (13.250 . it is direct to check by multiplying common 
factors and re-scaling the involved constants, particularly the angular wave numbers and the 
wave frequency, that all exact solutions presented in [^U] are special examples of combination 
solutions of the travelling wave solutions, defined in (13. 7p . (13. 8p . (I3.12p and (I3.13p . and some 
linear function solutions. 
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4 Concluding remarks 



A new systematical solution procedure to constructing exact solutions to nonlinear partial 
differential equations, both constant-coefficient and variable-coefficient, is proposed, based 
upon rational function transformations. Its key point is to search for rational solutions to 
variable-coefficient ordinary differential equations transformed from given partial differential 
equations. Together with an auto-Backlund transformation, an application of our method 
leads to various new travelling wave solutions to the 3-1-1 dimensional Jimbo-Miwa equation 

Our method allows us carry out the solution process of nonlinear wave equations more 
systematically and conveniently by computer algebra systems such as Maple, Mathemat- 
ica, MuPAD and Matlab, unifying the tanh-function method, the sech-function method, the 
homogeneous balance method, the extended tanh-function method, the sine-cosine method, 
the tanh-coth method, the Jacobi elliptic function method, the exp-function method, the 
F-expansion method, the mapping method and the extended F-expansion method. The pre- 
sented method also works for systems of nonlinear wave equations, and it can be applied to 
construction of special function solutions, like the extended F-expansion method [21], [22], the 
extended mapping method [39] and the improved extended F-expansion method [IQ] [H] [12] 

m- 

Besides travelling wave solutions and variable separated solutions, the 3-1-1 dimensional 
Jimbo-Miwa equation (11.11) has the following class of polynomial solutions: 

3 

'"1 = ^0 + ttix + + a^z + a^t + a^xz + a^xt + a^yz + -a^yt + a%zt, (4.1) 

where Oj, < i < 8, are arbitrary constants. These are all polynomial solutions found with 
Maple among the class of polynomial functions with individual degrees of the independent 
variables less than two. Letting u = Ui, the condition (13.161) becomes 

a2Vxx + aiVr^y = 0, (4.2) 

and thus, one can easily generate other solutions, for example, 

U2 = '"i|a2=o + f{x, t) + h{z, t), us = ui+ g{y, z) + h{z, t), (4.3) 

where ui is defined by ( 14.11) and /, g, h can be arbitrary functions. Taking /, g, h as polynomials 
in the indicated variables engender two other classes of polynomial solutions. 

All the presented solutions, including travelling wave solutions, variable separated solutions 
and polynomial solutions, show the remarkable richness of the solution space of the 3-1-1 
dimensional Jimbo-Miwa equation (11.11) . though the equation itself is expected to be non- 
integrable since it has solutions which are not single-valued functions in the neighborhood of 
their singularity surfaces [21]. A similar diversity situation of exact solutions can be found 
in solution spaces of typical nonlinear wave equations, for instance, the 1-1-1 dimensional 
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equations: the Korteweg-de Vries equation [IH HH], the Boussinesq equation [IHl HT], the 
nonhnear Schrodinger type equation [IH] and the Hirota-Satsuma coupled KdV equation 
and the 2 + 1 dimensional equations: the Kadomtsev-Petviashvili equation [50], the Davey- 
Stewartson equation [^11 E2] and the Boiti-Leon-Pempinelli dispersive long- wave system [H3] • 
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